A depression before a bump in the highest energy cosmic ray spectrum 
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We re-examine the interaction of ultra high energy nuclei with the microwave background radiation. 
We find that the giant dipole resonance leaves a new signature in the differential energy spectrum of 
iron sources located around 3 Mpc: A depression before the bump which is followed by the expected 
cutoff. 
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I. INTRODUCTION 

In 1966 Greisen, Zatsepin and Kuz'min noted that 
the microwave background radiation (MBR) makes the 
universe opaque to cosmic rays of sufficiently high energy, 
yielding a steep drop in the energy cosmic ray spectrum 
at approximately 5 x 10^^ eV (GZK cutoff). More re- 
cently, a fresh interest in the topic has been rekindled 
since several extensive air showers have been observed 
which imply the arrival of cosmic rays with energies above 
10^" eV. In particular, the Akeno Giant Air Shower Ar- 
ray (AGASA) experiment recorded an event with energy 
1.7 - 2.6 xl02° eV the Fly's Eye experiment re- 

ported the highest energy cosmic ray event ever detected 
on Earth, with an energy 2.3 - 4.1 xlO^o eV [|,§, both 
events being well above the GZK cutoff. Deepening the 
mystery, the identification of the primary particle in these 
showers is still uncertain. On the one hand, the Fly's Eye 
group claims that there is evidence of a transition from a 
spectrum dominated by heavy nuclei to one of a predom- 
inantly light composition ||^, while on the other hand, 
it has also been suggested that a medium mass nucleus 
also fits the shower profile of the highest energy Fly's Eye 
event 0. In addition, there is an unexpected energy gap 
before these events. Although heavy nuclei can be accel- 
erated to high terminal energies by "bottom up" mecha- 
nisms, one should note that, for energies above 100 EeV 
the range of the corresponding sources is limited to a few 
Mpc 1 11 1 . Sigl and co-workers § have analysed the struc- 
ture of the high energy end of the cosmic ray spectrum. 
They found that most "bottom up" models can be ruled 
out except for those involving a nearby source, which is 
consistent with data at the la level. Their argument for 
this is that a nearby source can account for the ultrahigh 
energy events but would also produce events in the ap- 
parent gap in data obtained to date. In this direction, 
Elbert and Sommers have suggested that the highest en- 
ergy event recorded by Fly's Eye, could have been ac- 
celerated in the neighborhood of M82, which is around 3 



Mpc away [^,^ . In relation to the aforementioned possi- 
bilities, we have re-examined the interaction of ultrahigh 
energy nuclei with the microwave background radiation 
and we have found a new feature in the ultrahigh energy 
cosmic ray spectrum from iron sources located around 3 
Mpc which forms the motivation for the present article. 



II. ENERGY ATTENUATION LENGTH OF 
ULTRAHIGH ENERGY NUCLEI 

The energy losses that extremely high energy nuclei 
suffer during their trip to the Earth are due to their in- 
teraction with the low energy photons of the MBR which 
they see as highly blue-shifted. The interaction with 
other radiation backgrounds (optical and infrared) can 
be safely neglected for nuclei with Lorentz factors above 
2 X 10^. Although the interactions of extremely high en- 
ergy nuclei with the relic photons lead to step-by-step en- 
ergy loss (which needs to be included in a transport equa- 
tion as a collision integral), in what follows we use the 
continuous energy loss approximation assuming straight 
line propagation which is reasonable for the energies and 
distances under consideration in this paper. The relevant 
mechanisms for energy losses are photodisintcgration and 
hadron photoproduction (which has a threshold energy 
of w 145 MeV, equivalent to a Lorentz factor of 10^^, 
above the range treated in this article) ||l^. 

The disintegration rate of a nucleus of mass A with 
the subsequent production of i nucleons is given by the 
expression |l3| . 
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dWr Wr(TAi{Wr) (1) 



where n{w) is the density of photons with energy w in the 
system of reference in which the microwave background 
is at 3K and Wr is the energy of the photons in the rest 
frame of the nucleus. As usual, F is the Lorentz factor 
and aAi is the cross section for the interaction. Using 
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FIG. 1. Energy of the surviving nuclei vs. propagation 
distance. It is also included the energy attenuation length of 
the surviving nucleon (dot line). 

the expressions for the cross section fitted by Puget et 
al. , it is possible to work out an analytical solution 
for the nuclear disintegration rates | [T5[ | . After summating 
them over all the possible channels for a given number of 
nucleons one obtains the effective nucleon loss rate. The 
effective ^^Fe nucleon loss rate obtained after carrying 
out these straightforward but rather lengthy steps can 
be parametrized by, 

R{T) ==3.25 X 10"'^r-"-^''3exp(-2.15 x 10^°/r) s'^ 

(2a) 

if r e [1. X 10^, 3.68 X 10^"], and 

R{T) = 1.59 X 10-12 p-0.0698 g-i (2b) 

if r e [3.68 X 10i",l. X 10"]. It is noteworthy that 
knowledge of the iron effective nucleons loss rate alone is 
enough to obtain the corresponding value of R for any 
other nuclei p^ ]. 

The emission of nucleons is isotropic in the rest frame 
of the nucleus, and so the averaged fractional energy loss 
results equal the fractional loss in mass number of the 
nucleus, viz., the Lorentz factor is conserved. The rela- 
tion which determines the attenuation length for energy 
is then, assuming an initial iron nucleus, 

i? = i?g e-^(r)*/56 (3) 

where Eg denotes the energy with which the nuclei were 
emitted from the source, and T = Eg/56. 

In Fig. 1 we have plotted the total energy of the heav- 
iest surviving fragment as a function of the distance for 
initial iron nuclei. Note that the values obtained here 




1e+20 1e+21 1e+22 



Eg 



FIG. 2. Relation between the injection energy of an iron 
nucleus and the final energy of the photodisintegrated nucleus 
for different values of the propagation distance (from grey to 
black 20 Mpc, 10 Mpc, 3.5 Mpc, 3 Mpc). 

are consistent with the ones obtained by Cronin using 
Monte Carlo simulation |ll]. One can see that nuclei 
with Lorentz factors above lO^" cannot survive for more 
than 10 Mpc 

In Fig. 2 the relation between the injection energy and 
the energy at a time t for different propagation distances 
is shown. The graph indicates that the final energy of the 
nucleus is not a monotonic function. It has a maximum at 
a critical energy and then decreases to a minimum before 
rising again as F rises, as was first pointed out by Puget 
et al. 0. The fact that the energy £^ is a multivaluated 
function of Eg leads to a pile-up in the energy spectrum. 
Moreover, this behaviour enhances a hidden feature of 
the energy spectrum for sources located beyond 2.6 Mpc: 
A depression that preceeds a bump that would make the 
events at the end of the spectrum (just before the cutoff) 
around 50% more probable than those in the depressed 
region. To illustrate this, let us discuss the evolution of 
the differential energy spectrum of nuclei. 

III. MODIFICATION OF THE COSMIC RAY 
SPECTRUM 

The photodisintegration process results in the pro- 
duction of nucleons of ultrahigh energies with the same 
Lorentz factor of the parent nucleus. As a consequence, 
the total number of particles is not conserved during 
propagation. However, the solution of the problem be- 
comes quite simple if we separately treat both the evolu- 
tion of the heaviest fragment and those fragments corre- 
sponding to nucleons emitted from the travelling nuclei. 
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The evolution of the differential spectrum of the surviv- 
ing fragments is governed by a balance equation that 
takes into account the conservation of the total number 
of particles in the spectrum. Using the formalism pre- 
sented by the authors in reference jl^, and considering 
the case of a single source located at to from the observer, 
with injection spectrum Q[Eg, t) — nE~^ S(t — to), the 
number of particles with energy E at time t is given by, 



kE-''+^ 

N{E,t)dE = — ^- dE, 

E 



(4) 



with Eg fixed by the constraint (||). 

Let us now consider the evolution of nucleons gener- 
ated by decays of nuclei during their propagation. For 
Lorentz factors less than 10"'^"'^ and distances less than 
100 Mpc the energy with which the secondary nucleons 
are produced is approximately equal to the energy with 
which they are detected here on Earth. The number of 
nucleons with energy E at time t can be approximated 
by the product of the number of nucleons generated per 
nucleus and the number of nuclei emitted. When the 
nucleons are emitted with energies above 100 EeV the 
losses by meson photoproduction start to become signif- 
icant. However, these nucleons come from heavy nuclei 
with Lorentz factors F > 10^^ which are completely dis- 
integrated in distances of less than 10 Mpc. Given that 
the mean free path of the nucleons is about A„ « 10 Mpc, 
it is reasonable to define a characteristic time Tj, given by 
the moment in which the number of nucleons is reduced 
to 1/e of its initial value Aq. In order to determine the 
modifications of the spectrum due to the losses which the 
nucleons suffer due to interactions with the relic photons, 
we assume that the iron nucleus emitted at t = to is a 
travelling source which at the end of a time has emit- 
ted the 56 nucleons together. In this way the injection 
spectrum of nucleons (F w 10^^) can be approximated 
by. 



0, 



(5) 



where Aq is the mass of the initial nucleus and the energy 
with which the nucleons is generated is given by Eq = 
Eg/Ao. 

The number of nucleons with energy E at time t is 
given by. 



n{E,t)dE 



kA, 



-7+2 



E, 



-7+1 



E 
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(6) 



and the relation between injection energy and the energy 
at time t remains fixed by the relation, A{t ~ t^) — 
Ei{B/E) + Ei{B/EG) = 0, Ei being the exponential 
integral, and A, B the parameters of the fractional energy 
loss of nucleons previously fitted by the authors [ p7[ . 

The modification factor 77, is defined as the ratio be- 
tween the modified spectrum and the unmodified one. In 
Fig. 3 we plot the modification factors for the case of 
sources of iron nuclei (propagation distance 20 Mpc) to- 
gether with the spectra of secondary nucleons. It is clear 
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FIG. 3. Modification factors for sources of iron nuclei at 20 
Mpc together with the spectra of secondary nucleons. 

that the spectrum of secondary nucleons around the pile- 
up is at least one order of magnitude less than the one 
of the surviving fragments. In Figures 4 and 5 we have 
plotted the modification factor for different propagation 
distances around 3 Mpc. They display a bump and a 
cutoff and, in addition, a depression before the bump. It 
is important to stress that the mechanism that produces 
the pile-up which can be seen in Figures 3, 4 and 5 is 
completely different to the one that produces the bump 
in the case of nucleons. 

In this last case, the photomeson production involves 
the creation of new particles that carry off energy yield- 
ing nucleons with energies ever closer to the photomeson 
production threshold. This mechanism, modulated by 
the fractional energy loss, is responsible for the bump in 
the spectrum. The cutoff is a consequence of the conser- 
vation of the number of particles together with the prop- 
erties of the injection spectrum (Jj^ E'g'^ dEg < 00). 

th 

In the case of nuclei, since the Lorentz factor is con- 
served, the surviving fragments see the photons of the 
thermal background always at the same energy. Then, 
despite the fact that nuclei injected with energies over the 
photodisintegration threshold lose energy by losing mass, 
they never reach the threshold. The observed pile-up in 
the modification factors is due solely to the multivalued 
nature of the energy at time t as a function of the injec- 
tion energy: Nuclei injected with different energies can 
arrive with the same energy but with different masses. 

It is clear that, except in the region of the pile-up, 
the modification factor 77 is less than unity, since r/ — 
{E / Eg)'^~'^ . This assertion seems to be in contradiction 
with the conservation of particle number. Actually, the 
conservation of the Lorentz factor implies. 



KEpdEgl = N{E,t)dE\, 



(7) 
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FIG. 4. Modification factor of single-source energy espectra 
for different values of propagation distance (from grey to black 
3 Mpc, 2.7 Mpc and 2.6 Mpc) assuming a differential power 
law injection spectrum with spectral index 7 = 2. 



FIG. 5. Same as Fig. 4 with spectral index 7 — 2.5. 
IV. CONCLUSIONS 



in accord with the conservation of the number of particles 
in the spectrum. Moreover, the condition (|^) completely 
determines the evolution of the energy spectrum of the 
surving fragments . Note that in order to compare the 
modified and unmodified spectra, with regard to conser- 
vation of particle number, one has to take into account 
that the corresponding energies are shifted. As follows 
from the conservation of the number of particles in 
the spectrum is given by, 

/ N{E,t)dE^ KEp dEg (8) 

with £^^1^ and £^7^1, ftic threshold energies for photodis- 
integration, and photopion production processes respec- 
tively. 

Let us now return to the analysis of Fig. 2 in relation 
to the depression in the spectrum. In the case of a nearby 
iron source, located around 3 Mpc, and for injection ener- 
gies below the multivalued region of the function E{Eg), 
E is clearly less than Eg and, as a consequence the de- 
pression in the modification factor is apparent. Then, 
despite the violence of the photodisintegration process 
via the giant dipole resonance, for nearby sources none 
of the injected nuclei are completely disintegrated yield- 
ing this unusual depression before the bump. For a flight 
distance of 3 Mpc, the composition of the arrival nu- 
clei changes from A = 50 (for F k, 10®) to A = 13 (for 
F K, 10^^). However, the most important variation takes 
place in the region of the bump, where A runs from 48 to 
13, being heavy nuclei of yl 33 the most abundant. For 
propagation distances greater than 10 Mpc one would ex- 
pect just nucleons to arrive for injection energies above 
9 X 10^° eV. In this case the function becomes multivalued 
below the photodisintegration threshold and then there 
is no depression at all. For an iron source located at 3.5 
Mpc, the depression in the spectrum is almost invisible 
(0(1%)), in good agreement with the results previously 
obtained by Elbert and Sommers using Monte Carlo sim- 
ulation g. 



We have studied the interaction of ultra high energy 
nuclei with the MBR. We have presented a parametriza- 
tion of the fractional energy loss for Lorentz factors up 
to 10^^ that allows us to analyse the evolution of the 
energy spectrum for different nuclei sources. When con- 
sidering an iron source located around 3 Mpc, the spec- 
trum exhibits a depression before a bump not previously 
reported. In the light of this finding it is tempting to 
speculate whether the apparent gap in the existing data 
is due to the relative weight of the depression and the 
bump if a source of iron nuclei is responsible for the end 
of the cosmic ray spectrum. This speculation, if true, 
reclaims "botton up" models as a possible scenario for 
the origin of the highest energy cosmic rays. The limited 
statistics in the observed data make it impossible to re- 
solve the question definitively at this time, and we are 
obliged to present this idea as a hypothesis to be tested 
by experiment. 

The existence of a cutoff or a gap which might be 
present in the observed spectrum is of fundamental in- 
terest in cosmic ray physics, allowing stringent tests of 
existing models. The future Pierre Auger Project jlj] 
should provide enough statistics for a final veredict on 
these open questions, and in particular on the ideas dis- 
cussed in this paper. 
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